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Abstract 

We revisit the WMAP dark matter constraints on Yukawa Unifi- 
cation in the presence of massive neutrinos. The large lepton mixing 
indicated by the data may modify the predictions for the bottom quark 
mass, enabling Yukawa unification also for large tan /3, and for positive 
\i that was previously disfavoured. As a result, the allowed parameter 
space for neutralino dark matter increases for positive particularly 
for areas with resonant enhancement of the neutralino relic density. 
On the contrary, a negative fi is not easily compatible with large lep- 
ton mixing and Dirac neutrino Yukawa couplings, and the WMAP 
allowed parameter space is in this case strongly constrained. 
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1 Introduction 



Reconciling the Cold Dark Matter (CDM) predictions of supersymmetric models with 
the stringent constraints from the Wilkinson Microwave Anisotropy Probe (WMAP) 
has been one of the challenges within the particle physics community in recent years. 
The amount of CDM deduced from the WMAP data is given by [1] 

n CDM h 2 = o.minis, (i) 

where h is the Hubble parameter in units of lOOkm/s/Mpc and flcDM = Pcdai/Pc the 
ratio of the matter density of cold dark matter, pcdm, over the critical density, p c , that 
leads to a flat Universe. This puts severe constraints on the type and parameter space of 
unified theories with Dark Matter candidates. One of the most promising frameworks 
in this respect is provided by supergravity models that conserve R-parity [2], and 
thus predict a stable lightest supersymmetric particle (LSP). This particle has to be 
neutral, and an obvious candidate is the neutralino [3], although gravitinos are also 
well-motivated [3] 0. 

If one imposes Yukawa unification, as expected from Grand Unified Theories (GUTs), 
the solutions become more predictive and additional constraints are imposed on the 
model parameters 0E]- The same holds for bounds from Flavour Changing Neutral 
Current (FCNC) processes, which have to be included in the analysis and strongly 
constrain the allowed supersymmetric parameter space [Hj- 

In addition to the above, the neutrino data of the past years provided evidence for 
the existence of neutrino oscillations and masses, pointing for the first time to physics 
beyond the Standard Model [9] . By now, it has been established that the atmospheric 
and solar mixing angles #23 and 9±2 are large and that the squared mass differences 
are Am 2 im ~ 2.5 x eV 2 and Am 2 oi ~ 8 x 10~ 5 eV 2 respectively [ID] . As ex- 
pected, however, the additional interactions required to generate neutrino masses also 
affect the energy dependence of the couplings of the Minimal Supersymmetric Standard 
Model (MSSM), and thus modify the Yukawa unification predictions. In fact, a first 
observation had been that the additional interactions of neutrinos to the tau would 
spoil bottom-tau Unification for small values of tan /3 where we are away from any 
fixed-point structure while the corrections to the bottom quark mass are small. Sub- 
sequently, however, it has been realized that the large lepton mixing could naturally 
restore unification, and even enable Unification for intermediate values of tan/5 that 
were previously disfavoured p21 02] • It is interesting to note that for the cosmologically 
favoured area, it is also possible to observe tau flavour violation at the LHC, in the 
framework of non- minimal supersymmetric Grand Unification [13] . 

In this work, we revisit the issues of Dark Matter and Yukawa Unification taking into 
account the effects of massive neutrinos, and large lepton mixing, as indicated by the 
data. As a first step, we extend previous results to large tan/3, finding significant 
effects on the allowed parameter space and on m&. In fact, it turns out that Yukawa 
Unification in the presence of neutrinos is also compatible with a positive /1 and large 
tan /?, unlike what happens if the effects of neutrinos are ignored [15], [161 E] • Passing to 
the relic density of neutralinos, we studied the consequences of the above, particularly 

*In fact, gravitinos can be dark matter even in theories with R-parity violation, if their lifetime is 
larger than the age of the universe [5] . 
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for the x — t coannihilation region and for resonances in the x~X annihilation channel, 
finding once more sizeable effects, which are further amplified when large lepton mixing 
is combined with quantum corrections above the GUT scale. 



2 Massive Neutrinos and Gauge and Yukawa Uni- 
fication 



The most straightforward extension of the Standard Model that can accommodate neu- 
trino masses is to include three very heavy right-handed neutrino states and assume 
that the smallness of masses arises from the See-Saw mechanism [18] . In this case, the 
predictions for mj and unification clearly get modified. In particular, radiative cor- 
rections from the neutrino Yukawa couplings have to be included for renormalization 
group runs from Mqut to Mx (scale of the heavy right-handed neutrinos) [19]. Below 
Mx, right-handed neutrinos decouple from the spectrum and an effective see-saw mech- 
anism is operative (with the neutrino mass operator running down to low energies). 
The relevant equations are given in [20] and are summarized in the Appendix. 

In addition, if the GUT scale lies significantly below a scale Mx, at which gravitational 
effects can no longer be neglected, the renormalization of couplings at scales between 
Mx and Mqut ma y induce additional effects to the running. The simplest such exam- 
ple is provided within the framework of the minimal supersymmetric SU(5) GUT, and 
the renormalization group equations (RGEs) in this case are given in [21J and are also 
summarized in the Appendix. These runs affect Yukawa unification via the supersym- 
metric corrections to the bottom quark mass; however, since in the leading-logarithmic 
approximation the induced modifications to soft masses are proportional to the Vqkm 
mixing [21], they are significantly suppressed. 

Nevertheless, it has been realized that the influence of the runs above the GUT scale 
on the Dark Matter abundance can be very sizeable [22] , due to changes in the relation 
between mf and m x , which is crucial in the coannihilation area. In dark matter calcu- 
lations, therefore, the possibility of such effects has to be also analyzed, and compared 
to the more standard scenarios. 

Our starting point will be the SU(5) superpotential at Mx which (omitting generation 
indices) is given by 



Here, T, F and S are the 10, 5* and 1 SU(5) superfields, respectively, and H and 
H are the 5 and 5* Higgs superfields. The couplings \ Ut d,N stand for the up-type 
quark, down-type quarks/charged lepton and Dirac neutrino Yukawa matrices. The 
symbol S stands for diagonal (the up- and down-type quark Yukawa matrices cannot be 
diagonalized simultaneously). Finally, Mx is the right-handed neutrino mass matrix. 
In addition, we work with the soft SUSY-breaking Lagrangian 



W x = T T X S U T H + T T X d F H + F T X 5 N S H + S T M N S 



(2) 




f f m 2 10 f + F f ml F + S f m\ S + 
m 2 h h j h + m\h ] h + M 5 \ 5L \ 5L + h.c. + 

{f T A u f h + f T A d F h + F T A N S h + h.c. } 



(3) 
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where m 10t51 are the 10, 5* and 1 SU(5) superfields masses, respectively, and m h ^ 
are the 5 and 5* Higgs superfields masses. A U ^ N are the up-type quark, down-type 
quark/charged lepton and Dirac neutrino trilinear terms. Finally, X^l is the SU(5) 
gaugino and M 5 its soft Majorana mass. In the mSUGRA scenario, universal boundary 
conditions are assumed at the gravitational scale 

m io = m\ = m\ = m h = m- h = to 2 ,, (4) 

A f = a \ f , f = {u,d,v} (5) 

The physical values of the masses will then be obtained by integrating the renormal- 
ization group equations from Mx down to low energies. In the first part of the runs 
we work with the RGEs of SU(5) from the high energy scale Mx down to Mqut- Sub- 
sequently, we run the MSSM RGEs, supplemented with right-handed neutrinos, from 
M GUT to the right-handed neutrino scale, M N . At this scale, right-handed neutrinos 
decouple from the spectrum, leaving us with the MSSM RGEs and the see-saw operator 
(which is also renormalized down to the low energy scale). The effects of neutrinos on 
Yukawa unification crucially depend on the magnitude of the dominant neutrino Dirac 
Yukawa coupling, A at, and consequently the magnitude of (which determines Xn 
through the See-Saw conditions). Obviously, for low M N and Ajv, the effects are less 
significant. 

Large neutrino Yukawa couplings affect unification by increasing the predicted value of 
m,b(Mz)- This can be understood for small tan/5 by simple, semi-analytic expressions 
[T2l [T3] since only the top and the Dirac- type neutrino Yukawa couplings (X t and A^) 
may be large at the GUT scale. In this case, the RGEs take simple form 

167r 2 ^A t = (6A 2 +X 2 N - G v ) A t , lQ ^J t X N = (4A^ + 3A 2 - G N ) X N 

lQ* 2 j t h = (A f 2 - G D ) X b , 16vr 2 ^A T = {X% - G E ) X T (6) 

Here, A a , a = U, D, E, N, represent the 3x3 Yukawa matrices for the up and down 
quarks, charged lepton and Dirac neutrinos, and G a = Yli=i c a9i(t) 2 are functions of 
the gauge couplings with the coefficients c^'s as in [llj. Denoting by Xb , X T0 the b and 
t couplings at the unification scale, it is found that 

A 6 (tAr)=p6— K(t N ), P=tV ( 7 ) 

IE A T0 4at 

where 7«(t) and depend only on gauge and Yukawa couplings. Let us then assume 
successful b — t unification at M GUT , with A ro = Xb . In the absence of right-handed 
neutrinos ^ N = 1, thus p — 1. In the presence of them, however, A ro = Af, at the 
GUT scale implies that p ^ 1 (since £at < I)- To restore p to unity, a deviation 
from b — t unification would seem to be required. However, large lepton mixing is 
reconciled with Yukawa unification, by making the simple observation that the b — r 
equality at the GUT scale refers to the (3, 3) entries of the charged lepton and down 
quark mass matrices. It is then possible to assume mass textures, such that, after 
the diagonalization at the GUT scale, the (m^ a9 ) 33 and (m^ a9 ) 33 entries are no-longer 
equal [12]. The simplest example the one with symmetric mass matrices: 

M^A( y M, MWa( X 1 l) (8) 
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Here, A may be identified with mb^McuT)- While the textures ensure equality of 
the (3,3) elements of the down and charged lepton mass matrices, the eigenvalue of 
the charged lepton mass matrix is not 1, but 1 + x 2 , thus implying that A& ^ X T 
after diagonalization. Within this framework, the issue of Yukawa Unification in the 
presence of large lepton mixing has been analyzed in [T2"l [T3] for small and moderate 
values of tan [3. 

In our current work, we extend these results to all values of tan/5 and both signs 
of fi, taking appropriately into account the large supersymmetric corrections to mj, 
[231 1211 125] . and using up-to-date experimental bounds. Potentially large lepton mixing 
effects will be taken into account by imposing the GUT condition A r = A{>(1 + 5) where 
S is determined by requiring that m^iMz) is in the allowed experimental range. The 
third generation Dirac neutrino Yukawa coupling A^ is determined through the see- 
saw mechanism by assuming m Vi = 0.05 eV and a heavy Majorana neutrino scale 
Mn = 3 x 10 14 GeV (a value that ensures that Xn stays within the perturbative 
regime). On the other hand, the behaviour associated with the absence of a Dirac 
neutrino coupling should be recovered by appropriately lowering M N and thus A at. 
The resulting effects to the allowed parameter space for unification are significant, 
and give interesting information on the magnitude and origin of lepton mixing. We 
then proceed to discuss in detail the cosmological implications in the framework of the 
WMAP data. 

3 Yukawa Unification in SUSY Models 

While gauge unification is considered as one of the most attractive features of Su- 
persymmetry, the relations among Yukawa couplings derived by embedding SU(3) x 
SU(2) x U(l) in a larger gauge group are less clear. The charged fermion mass hierar- 
chies indicate that any mass correlations induced by a unified gauge symmetry will be 
more explicitly manifest in the relation between the Yukawa couplings of the third gen- 
eration, X t , A;,, A T . How this works depends on the theoretical framework that defines 
the initial conditions. In 5*0(10), for instance, all fermions are in the same represen- 
tation, thus, in the simplest realizations, neutrino couplings will be unified with the 
rest. In SU(5), the field structure is (Q,u c ,e c )i G 10 of SU(5) and (L,d c )i e 5, only 
implying b — r unification; neutrinos are singlets and thus there is a freedom of choice 
(although the fact that they couple to the same Higgs as the up-type quarks could be 
a motivation for some link with the top coupling). Clearly, the more restrictive the 
unification schemes, the stronger the correlations between quarks and leptons. 

In supersymmetric models, unification turns out to be very sensitive to the model 
parameters. Among others, the compatibility of b — r unification with the observed b 
and r masses has a sensitive dependence on the sign of the Higgs mixing parameter, 
/i, and on the details of the superparticle spectrum [26]. Moreover, the universality 
condition on the soft terms can be removed as in [7], resulting in an enhancement of 
the allowed parameter space. In fact, it turns out that full third generation Yukawa 
Unification X t = X b = X T in the CMSSM fails to accurately predict the experimental 
values of the third generation fermion masses, unless one goes to a large tan (3 regime 
with a heavy spectrum and a dark matter abundance that would tend to overdose the 
universe. 
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To correctly obtain pole masses within this framework, the standard model and su- 
persymmetric threshold corrections have to be included; for the bottom quark, these 
corrections result to a Amj that can be very large, particularly for large values of tan f3. 
In mSUGRA, in the absence of phases, Amj has the sign of /x [231 |2"4"1 125] , which is 
required in order to obtain a 6-quark mass in the allowed range in models with b — r 
unification. However, this also results to a positive supersymmetric contribution to 
BR(6 — > 57), which can be compatible with the bounds only for a heavy sparticle 
spectrum. 

The allowed parameter space is nevertheless extremely constrained from the bounds 
on Flavour Changing Neutral Currents, and the new bounds on b — > 37 are crucial 
for the whole discussion and comparisons with the SM prediction [25] . The — 2 
constraints from Brookhaven National Laboratory (BbNL) [29] are also relevant since 
the supersymmetric contribution to the muon anomalous magnetic moment takes the 
sign of fj, in mSUGRA [30]. There are also some uncertainties 0(a 2 ) in the Standard 
Model prediction, due to the hadronic vacuum polarization correction [3T]- For a 
heavy sparticle spectrum the SUSY contribution to Aa M is small with respect to the 
Standard Model one. However, we do not take into account that constraint derived from 
g M — 2, since there is no significant deviation from the SM predictions if the hadronic 
contribution is calculated using r-decay data (the current experimental value deviates 
by 3.4 a from the SM prediction if the hadronic vacuum polarization is determined 
using e + e~ annihilation data). 

Let us also summarize a few facts on the possible range of the mass of the bottom 
quark: The 2-a range for the MS bottom running mass, m b (m b ), is from 4.1-4.4 GeV 
(with corresponding pole masses from 4.7 to 5 GeV). We also know that a s (Mz) = 
0.1172 ± 0.002, and the central value of a s corresponds to m b {Mz) from 2.82 to 3.06 
GeV, while the value m b {m b ) = 4.25 GeV is mapped to m b {Mz) = 2.92 GeV. The 
allowed strip for m b (M z ) moves to 2.74 GeV < m b (M z ) < 3.014 GeV for a™ ax and to 
2.862 GeV < m b (M z ) < 3.114 GeV for a™ n . 

In the left panel of Fig. [TJ we summarize the predictions for m b in mSUGRA, where 
all CP phases are either zero or 7r (defining the sign of /i). In order to discuss the 
dependence of m b [Mz) on tan/3, we consider the following set of soft parameters: 
Mi/a = 800 GeV, A = 0, m = 600 GeV. We study this set for both \i > and 
/i < 0, setting a s (M z ) = 0.1172. We also include a reference line without the SUSY 
corrections to the bottom mass (double-dot-dash line, for Am;, = 0). The figure 
exhibits the known fact that in the absence of phases or large trilinear terms, Am b 
is positive for /1 positive, and therefore the theoretical prediction for the b quark pole 
mass is too high to be reconciled with b — r unification. On the other hand, for < 0, 
Am b is negative and the theoretical prediction for the 6-quark mass can lie within the 
experimental range for values of tan f3 between roughly 30 and 40; clearly, for a large 
tan/3 it is mandatory to take into account the large supersymmetric corrections to 
m b [231 [2H I2S]- The figure also illustrates that, for ji > and after taken properly into 
account supersymmetric corrections, m b scales very slowly with tan f3. In other words, 
the renormalization flow for the bottom mass is attracted by some fixed value regardless 
of the initial conditions [2HE2] (which as we see in the figure turns out to be too large). 
This quasi- fixed point, however, is not generated purely by the renormalization flow 
given by the RGE, and only appears after including supersymmetric corrections (the 
curve for Am b = is not flat in Fig. [T]). Furthermore, for /i < 0, there is no fixed point 
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at all. 
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Figure 1: The value of m^Mz) versus tan/3 assuming A& = A T at the high scale in the 
absence (left plot) and presence (right plot) of massive neutrinos and for both signs of 
H for the following set of parameters: Mi/ 2 = 800 GeV , A = GeV, m = 600 GeV. 
The experimental range of m b (horizontal lines) shown in the left pannel is computed 
for the central value of a s (Mz). In the right panel, the solid lines are obtained within 
the MSSM, while the dot-dash lines include Dirac Yukawa couplings up to the scale 
Mn = 3 x 10 14 GeV. Here we take the lowest (blue thin) and highest (black thick) 
experimental values of a s (Mz). 



In the right panel of Fig. [T]we repeat the analysis in the presence of massive neutrinos, 
keeping only the third generation couplings (and ignoring lepton mixing effects) from 
the Mqut to the scale of the right-handed neutrino masses, Mat, and evolve the light 
neutrino mass operator from this scale down to Mz- A large value of the Dirac- 
type neutrino Yukawa coupling, Ajy, at the GUT scale may arise naturally within 
the framework of Grand Unification, and its value is determined, through the see- 
saw mechanism, by demanding a third generation low energy neutrino mass of m U3 = 
0.05 eV and evaluating the respective heavy neutrino Majorana mass from the see-saw 
conditions. The predictions for rrib(Mz) using the lower and upper bounds of the 2-cr 
experimental range of a s and the corresponding range for m^Mz) after the evolution 
of the bounds on m^m^) are shown for M^ = 3 x 10 14 GeV. We use m t = 172.6 GeV 
[3"3] in all computations. 

We observe that for \i > the prediction for m^Mz) is always very large, despite its 
dependence on the soft terms through Amj. For /i < 0, there is a window of values of 
tan j3 compatible with asymptotic b — t Yukawa unification. For the values of the soft 
terms considered in Figfj], the allowed range of tan f3 moves from 27 — 44 to 30 — 45 
when we introduce the effect of see-saw neutrinos. Let us stress that b — r Yukawa 
unification is (is not) compatible with m&, for fi < (> 0), regardless on whether or 
not massive neutrinos are included, when lepton mixing effects are ignored. 
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4 Yukawa Unification and for large lepton mix- 
ing 



The results are significantly modified once we consider the effects of lepton mixing in 
the diagonalization and running of couplings from high to low energies. In order to show 
this, we focus on b — r unification within the framework of SU(5) gauge unification and 
flavour symmetries that provide consistent patterns for mass and mixing hierarchies, 
and naturally reconcile a small Vqkm mixing with a large charged lepton one. Taking 
into account the particle content of SU(5) representations (with symmetric up-type 
mass matrices, and down-type mass matrices that are transpose to the ones for charged 
leptons), one finds that 

•M u oc(J E [ ^, Ml* A (I J), MWA^l *) (9) 
which, after diagonalization, lead to the relation 

ml m!l „ o / i .-) ,.2 



1 + x 2 1 — x 1 



1 - 2a; 2 +0 (6 2 ) (10) 



where 5 parametrizes the flavour mixing in the (2,3) charged lepton sector (any addi- 
tional mixing required to match the data would then arise from the neutrino sector). 

In the left panel of Fig[2] we show the change of m& as a function of tan /3, when the 
effects from large lepton mixing are correctly considered. Comparing with the previous 
plots, we see how solutions with positive \x are now viable, for the whole range of tan f3. 
The appropriate size of the parameter 5 in each case can be determined by imposing 
the relation A T = A&(1 + 5) at Mqut and investigating the values that are required in 
order to obtain a correct prediction for m&(Mz). This is shown in the right panel of 
Fig 121 where we demand a value of m^iMz) at the center of its experimental range, for 
the central value of a s . We checked (although not shown) that the effect from runs 
above M GUT is very small, as expected. We then observe the following: 

• Solutions with a positive /i are now enabled, for the whole range of tan j3. 

• Negative \x is also allowed in the whole range of tan/3 and requires a smaller 
mixing parameter 5 than the \x > case. 



5 Yukawa unification and Dark Matter constraints 

In mSUGRA (or the CMSSM) for choices of the soft terms below the TeV scale, the 
LSP is mainly Bino-like and the prediction for Q x h 2 is typically too large for models 
that satisfy the experimental constraints on SUSY. These constraints exclude models 
with relatively small values of mo and My2 in which neutralino annihilates mainly 
through sfermion exchange. The values of WMAP can be obtained mainly in two 
regions: 

• \ — f coannihilation region, which occurs for m x ~ m f . 
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Figure 2: In the left panel, we show as a function of tan (3 for /i > 0, without (solid 
lines) and with (dot dashed) massive neutrinos, in the latter case assuming m Uz = 0.05 
eV and Mjq = 3 x 10 14 GeV. In the right panel, we show the required values of 5 
compatible with b — r unification at the central experimental value of for both signs 
of n- 



• Resonances in the x~ X annihilation channel, which occur for m,A ~ 2m x . 

Since the above areas are tuned [31], they will inevitably be sensitive to the changes 
induced by GUT unification and sizeable mixing in the charged lepton sector. In this 
respect, it is illustrative to first investigate the impact of the parameter 5 on the value 
of m A . We show this in Fig. [3J, where we present the variation of rriA with tan/3 for 
m = 600 GeV, M 1/2 = 800 GeV and A = 0. In the upper lines the value of 0) 
is fixed by demanding m^iMz) = 2.921 GeV while in the lower ones 5 = with no 
restrictions on the m&(Mz) prediction. We observe that the upper lines do not meet 
the resonance condition, while for the lower lines this condition is achieved for values 
of tan/3 in the range 40 — 50. We also see that, assuming soft term universality at a 
scale Mx > Mqut, the runs beyond Mqut change only moderately the prediction for 
m A . 

Our next step is a global study of the supersymmetric parameter space by assum- 
ing universal soft terms m , M 1 / 2 and A at some scale M x > M GUT , and obtain 
the mass spectrum by integrating the appropriate RGEs at each energy range, as de- 
scribed above. The physical observables are computed using the code provided by 
micromegas [35]. The 2-cx range for BR(6 — > 57) is constructed including an intrinsic 
0.15 ■ 10~ 4 MSSM correction as in Ref. [36], leading to 

2.15 • 10~ 4 < 5^0 -> 57) <4.9-10~ 4 (11) 

Without including the MSSM intrinsic error, the above range becomes 

2.8 ■ 10~ 4 < BR(b -f S7) < 4.4 ■ 10~ 4 (12) 
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Figure 3: Values of the pseudoscalar higgs mass tua versus tan(3 for the same setting 
of parameters as in FigUl In the upper lines the parameter 6 varies so that rrib(Mz) = 
2.921 GeV, while in the lower lines 5 = 0. The solid lines correspond to Mx = 
2 • 10 17 GeV and the dashed lines to Mx = Mqut- 



We take a conservative bound for the mass of the lightest CP-even Higgs of = 
114 GeV and also consider an uncertainty of ~ 3 GeV [37] to its theoretical computa- 
tion. 




M 1/2 (GeV) M 1/2 (GeV) 

Figure 4: WMAP allowed area (green-shaded) for the case of tern (3 = 45, /i > 0, A = 0, 
when a s (Mz) and rrib(Mz) are set at their central values. The solid black lines indicates 
the b -»• S7 constraints [thick BR(b -> 57) = 2.15 • 10~ 4 , thin BR(b -> 57) = 2.8 • 10" 4 /, 
and the dash blue line the Higgs mass bound [thinm h = 114 GeV, thickm h = 111 GeV]. 
On the left graph, M x = M GUT while on the right panel, M x = 2 ■ 10 17 GeV. 
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In Fig. H]we present the WMAP favored area for tan/? = 45, including lepton mixing 
effects and considering the most tolerant bounds on BR(6 — > 57) and m/ 4 (stricter 
bounds on these constraints are also displayed). In this example, we keep m^iMz) = 
2.920 GeV which corresponds to the evolution of the experimental value of m^ obtained 
for the central value a s (M z ) = 0.1172. When M x = M GUT we find the WMAP favored 
region to be on the familiar CMSSM x~ T coannihilation area. However, as in [22], we 
find that the runs corresponding to Mx > Mqvt have a big impact on the neutralino 
relic density and Fig. H] shows clearly this effect. The large values of the gauge unified 
coupling asu(5) tend to increase the values of m f in a way that, even if we start 
with mo = at Mx the model predicts > m x . This implies that the allowed 
parameter space with a neutralino LSP is significantly enhanced (green area), while 
the coannihilation condition becomes harder to achieve; this effect is more visible for 
large tan (3 and is understood since the lightest stau has a mass m fl ~ m fRR + m fLR , 
where 

m f R R — (1 - Pp) m 2 + 0.3Mjy 2 , m fLR ~ -m ry utan/3 (13) 

P/3 being a positive coefficient dependent on tan/? (in our case pp < 1) [22]. The 
increase of m ?1 is due to the GUT runnings (~ M^ 2 ). The picture for tan[3 = 35 
and A = m has been discussed in [13], where it was shown that the WMAP allowed 
region can be compatible with observable flavour violation at the LHC. 

We stress again that b — r Yukawa unification in the context of the CMSSM requires 
5 > 0. An exhaustive study of the parameter space for the case of Aq = and keeping 
the prediction of m^Mz) = 2.92 GeV is presented in Figs. [5] and El The left panel 
of Fig. [51 shows that the phenomenological and cosmological constraints are fulfilled 
for tan/3 > 31 and Mx/2 > 330 GeV respectively. On the other hand, the right panel 
of Fig. [5] indicates that sizeable values of the mixing parameter 5 are needed in order 
to maintain rrib(Mz) in the experimental range and that the required mixing increases 
with tan/3. Fig. [6] shows the (mo, M 1 / 2 ) plane for different values of tan/3, taking 
all constraints into account; it can be seen that the resonant effects start becoming 
important at tan/3 > 45. 

The impact of varying m;, and a s within their experimental range can be seen in Fig. [7] 
which indicates that, for tan/3 below 40 there are no significant changes with the 
variation of on the WMAP area that lies on the coannihilation region. For larger 
values of tan/3, however, small changes of the bottom Yukawa coupling due to the 
modification of rrib have a significant impact on resonant annihilation, as is manifest 
by the locations of the bands for tan f3 = 45 and 50 in the two lower plots. 

The case /i < is not yet ruled out by the (g — 2)^ data, since as we mentioned before 
by using r data, a small negative discrepancy of the experimental measurement as 
compared to the SM prediction can be accommodated. However, the upper limit on 
the BR(6 — > S7) imposes a severe constraint on the SUSY parameter space since the 
supersymmetric contribution adds to the SM prediction. 

For p < and Mx > Mqut-, we find that the WMAP areas allowed due to \ — f 
coannihilations at low values of tan/3 are excluded by the b — >• 57 higher bound. 
However, at larger values of tan /3, the areas with resonant annihilations are marginally 
allowed, as we can see in Fig. El 

As already underlined, these results are very sensitive to Xx- Once we decrease it 
by lowering Mx (so that, from the see-saw condition, m U3 remains 0.05 eV) we allow 
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M 1/2 (GeV) tan (3 

Figure 5: The left panel provides a plot equivalent to Fig. 4 of |Hj/ for the purpose of 
comparisons; the upper part corresponds to areas with resonant Higgs channels. In the 
right panel, we show the values of 5 corresponding to the previous plot; these inevitably 
implies 5^0. 




Figure 6: WMAP allowed areas for several values of tan (3. We consider m^iMz) — 
2.92 GeV, A = GeV and fi > 0. 

larger regions of the parameter space for negative \x as well. In both cases, however, 
the values of the parameter 5 have to be very small (5 < 0.05); therefore, negative \x is 
not compatible with large charged lepton Yukawa mixing. 

6 Conclusions 

We revisited the WMAP dark matter constraints on Yukawa Unification in the presence 
of massive neutrinos. Large lepton mixing, as indicated by the data, modifies the 
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Figure 7: Same plot as Fig® but taking the extreme values of both a s (M z ) and 
mb(Mz). On the leftpanel, a s (M z ) = 0.121 andmb(Mz) = 2.74 GeV while on the right 
panel a s (M z ) = 0.113 and m b (M z ) = 3.114 GeV. In both cases M x = 2 ■ 10 17 GeV. 




M 1/2 (GeV) M 1/2 (GeV) 

Figure 8: (m , M 1 / 2 ) plane for tan(3 = 40, A = 0, and /i < . On the left panel 
we take the right-handed neutrino scale = 6 x 10 14 GeV, while on the right panel 
M N = 10 13 GeV. Here, m b (M z ) and a s (M z ) are set to their central values. The lines 
follow the same notation as Fig. [^} 

predictions for the bottom quark mass, and enables Yukawa unification also for large 
tan/? and for positive values of fi, which were previously disfavoured. The larger the 
Dirac neutrino Yukawa couplings, the larger the effects. A direct outcome is that the 
allowed parameter space for neutralino dark matter also increases, particularly when 
the effects of large lepton mixing are combined with runs above Mqut- Summarising, 
we find the following: 
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• b — r unification is only allowed for \i > in the presence of large charged lepton 
mixing, which is motivated by the experimental data of the recent years. 

• For n < we still find values compatible with an exact b — r unification at 
Mx- However, for large Dirac neutrino Yukawa couplings and mixing arising 
dominantly from the charged lepton sector in the basis where the down-quark 
mass matrix is diagonal, the space of parameters compatible with WMAP is i) 
only marginally compatible with the upper bound on b — > sj and ii) entirely 
excluded on the grounds of the (g — 2) M observations, if we use e~e + data to 
estimate the SM vacuum polarization contribution. 

Interestingly enough, it turns out that the cosmologically favoured parameter space also 
implies lepton flavour violating rates that are very close to the current experimental 
bounds [38] . Finally, additional interesting effects may arise in the case of non-universal 
soft terms; these are also addressed in detail in [38] . 
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Appendix 

In this appendix we summarize the RGEs that are most relevant for the purposes 
of the work addressed in this paper. For runs above the GUT scale the equations 
involving the Yukawa couplings and the soft mass terms corresponding to the 10 and 
5 representations of £77(5), for the 3 rd generation, take the form 



167T 



2 d\]\[ 



16tt' 



16tt 



dt 
2 d\t 



dt 



.^-gl + 7X 2 N + 3X 2 t +4X 2 b 



84 



g 2 5 + 10X 2 d + 3X 2 t +X 2 N 



96 



9l + 9A? + 4A 2 + A 2 



N 



X 



N 



x d 



(14) 
(15) 
(16) 



16tt 2 



dm 



10 



144 



dt 



^ 5 2 M 2 +(12A 2 + 4A 2 ) 



m 



10 



+4 [(ml + ml) A 2 + A 2 d ] + 6 (A 2 m\ + A 2 ) 



(17) 



16vr 2 ^! 
dt 



96 2 lJ 



+ A 



I) 



2 (4A 

+8[(m 2 10 + m 2 h ) X 2 d + A 2 d ] 



m r 



;is) 



2 (A 



N m h 



+ A 



N m l 



+ A 



14 



For runs from Mqut to Mat, the equations for the Yukawa matrices are: 



16tt z 



dX 



N 



dt 



(4X% + 3A? + XI) 



X 



N 



(19) 



16tT 



dX T 
~dT 



: gj + 3g 2 2 ] h 



(4A^ + 3A^) 



A T 



16-7T' 



dXt 
dt 
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9l + 3^ 2 2 + y^ 3 2 ) h - (6A? + Ag) + A 



A, 



(20) 



(21) 



Since the neutrino has no coupling to the bottom quark, the Yukawa matrix corre- 
sponding to the latter remains unchanged with respect to the MSSM case. 

In section 2 semi-analytical expressions for the small tan f3 regime were given. In that 
case, only the top and the Dirac-type neutrino Yukawa couplings can be large at the 
GUT scale. 
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